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Abstract. We investigate the electronic and superconducting properties of a 
ncgativc-l/ Hubbard model. For this purpose we evaluate a recently introduced 
variational theory based on Gutzwiller-correlated BCS wave functions. We find 
significant differences between our approach and standard BCS theory, especially 
for the superconducting gap. For small values of \U\, we derive analytical express- 
ions for the order parameter and the superconducting gap which we compare to 
exact results from perturbation theory. 
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1. Introduction 

Gutzwiller-correlated wave functions pQ were originally proposed for a variational 
examination of the one-band Hubbard model j^j 

H = J2 ^'fitrCjo + U Y1 ^T^i = #kin + H V (1) 

i¥=j 0-=T,l i 

with positive Coulomb repulsion U. Here, i and j are sites of the respective lattice 
under investigation and cf a creates (annihilates) an electron with spin a =T,|- 

As usual, hi,, = cf a Ci a counts the number of cr-electrons on site i. 

Gutzwiller's primary intention was to study ferromagnetism in such a system. 
Unfortunately, the evaluation of expectation values for Gutzwiller's variational 
wave function still poses a generally unsolvable many-particle problem. Gutzwiller 
introduced an approximate evaluation scheme, the so called Gutzwiller approximation, 
which was based on phenomenological counting arguments; for a mathematically 
sound introduction to this technique, see More recent research on Gutzwiller 
wave functions focused on the development of better controlled evaluation schemes. 
Main progress was made in the limit of infinite spatial dimensions 0J |S] where the 
Gutzwiller approximation turned out to be exact in some cases UJ. 

Despite its general importance as the most simplistic model for correlated electron 
systems, the one-band Hubbard model is only of limited use for the investigation of 
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real materials. Therefore, a second aim was the evaluation of more general Gutzwiller 
wave functions for multi-band models. Such a theory has first been derived in [7] 
and was later used in numerical studies on Nickel jS] EJ- The results of these nu- 
merical investigations are in excellent agreement with experiments and pose a major 
improvement over band-structure calculations based on the local-spin-density approxi- 
mation to density-functional theory. 

Since the advent of high-temperature superconductivity about 20 years ago, 
Gutzwiller wave functions have been used for the investigation of systems with alleged 
attractive interactions, most notably the two-dimensional i-J-model 111! IT2*| . For 
the evaluation of those Gutzwillcr-correlated BCS wave functions approximations 
have been used which were based on counting arguments similar to the original 
Gutzwiller approximation. In [Hj we developed a general theory for the evaluation of 
superconducting multi-band Gutzwiller wave functions in infinite spatial dimensions. 
An application of these general results to real materials requires significant numerical 
efforts. In this report we present first results on the simplest correlated electron 
system that exhibits superconductivity, the attractive Hubbard model. The same 
system has been studied using Gutzwiller-correlated wave functions in pi], based 
on an evaluation scheme first derived in For the attractive one-band Hubbard 

model with a Gaussian density of states, the numerical results of both methods seem 
to agree. Unfortunately, we cannot provide a rigorous proof for the equivalence of 
both approaches. Since the evaluation scheme used in ^] cannot be generalized 
easily to more realistic model systems, we do not further pursue their approach. 

We structure our article as follows. In section we introduce the Gutzwiller- 
correlated superconducting wave functions for the negative- U Hubbard model and 
derive analytical expressions for the variational ground-state energy. The minimisation 
of this energy is discussed in section[3] In particular, we derive an effective one-particle 
Hamiltonian which provides a link to ARPES experiments. In section 0] we focus on 
the half-filled Hubbard model. We discuss the small-] U\ limit analytically and provide 
numerical results for a system with a semi-elliptic density of states. A summary in 
section [S] closes our presentation. 



2. Superconducting Gutzwiller wave functions 

We consider the Gutzwiller wave function 

|¥o> = Ip*l*°>- W 

i 

Here, |^o) is a normalized quasi-particle vacuum and the local correlator 

Pi = P = Y l Xi,u\I)ii(I'\ (3) 

induces transitions at the lattice site i between the four atomic states \J), i.e., the 
empty state |0), the doubly occupied state \d) and the two single electron states 
\ a ) = I T); I I)- We denote expectation values with respect to \^g) and |*o) as (■ • -)g 
and (• • -)o, respectively. In general, there are 16 variational parameters Xij> . When we 
assume that there is no magnetic or charge order in the system, it is sufficient to work 
with only four independent and real parameters: (i) Ag = (ii) Ai = Af.f = Aj_. j_ , 
(iii) Xd = Xd.d, (iv) Xb = A<2,0 = X® t d- It can be shown that a more general Ansatz 
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with complex parameters Aj i> does not lead to a variational improvement. As shown 
in 9 the four parameters are not independent. Instead, they obey the constraints 



i =(Pr 



no = {cfP 2 c ia ) , (4) 
Ac = (c+P 2 c+) . 

Here, no and Ao are the elements of the uncorrelated local density matrix 

"0 = (c+ c ia } , A = (c+ T c+) (5) 
which we assume to be independent of the site index i and the spin index a. The 
evaluation of equations leads to 

1 = (A 2 , + A|)m 0)O + 2A?mi, + (A 2 , + \%)d + 2A O (A + X d )X B , 
no = (A| + \ 2 B )mi,o + \ 2 d , (6) 
A = A 2 A - (A + A d )A B mi !0 , 
with the expectation values 

d = (h^nii)o = n 2 a + Aq , 

mi, = - fiii))a = n Q - d , (7) 

m0, o = ((1 - rV)(l - «a))o = 1 - 2mi,o - • 
We use these equations in order to express the three parameters A , Ai, A B in terms 
of the fourth parameter A^, 

A x = 1 + x(d — n ) , 

Xl = X 2 d - x(l - 2n ) , (8) 
X%=l-X d + a;m , o . 
Here, we introduced the abbreviation 



B I I AA 2 C 

x = - 

with 



X =2A 2 V 1 -^- 1 (9) 



A = A 2 + (1 - 2n o )m 0iO , 

B = — 4A 2 ,A 2 1 m , o + 2A(1 ~ A 2 )(l ~ 2n ) , (10) 

C = (1 - A 2 ) ((1 - 2n ) 2 (l - X 2 ) - 4A 2 A 2 ) ■ 
The expectation value of the Hamiltonian Q for the Gutzwiller wave function (J2J 
can be evaluated in the limit of infinite spatial dimensions UJ. For a hopping term 
this evaluation yields 

+ 99"((ctc+ ff ) + {^v)o) (11) 

where 

<7 = q(X d ,n , A ) = Ai(A (l - n ) + A d n ) + 2A B A Ai , (12) 
q = q(X d , n , A ) = A Ai(A d - A + A B )(1 - 2n ) (13) 

are renormalization factors for normal and anomalous hopping. The expectation value 

of the Hamiltonian then reads 

E™ = (H) G = Q(* |ffkin|*o) + Ud = QE° in + Ud , (14) 

with 

Q = q 2 -q 2 , d = X 2 d d + X 2 B m^ + 2X B X d A Q . (15) 
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3. Minimization procedure 

The variational energy has to be minimized with respect to Ad and |\&o) whereby the 
equations 101 and the normalization of |^o) need to be obeyed. For this purpose we 
introduce Lagrange parameters r/ n , i] s and E SP . Furthermore, we keep the average 
number of particles 

n=-T = n(\ d , n , A ) = (2A?mi, + 2d) (16) 
fixed by means of a Lagrange parameter u. The minimisation problem then becomes 
££ ar = Minimum \E val - E SP ((* |*o> - 1) 

- ^(A„ - (*o|S+ T 6+ |* > + cc.) (17) 

- Vn ^2(n Q - {^ \ct,Ji,a\^o)) + (*L(n - n(X d , n , A )) 

i,(7 

The minimization with respect to |\&o) can be carried out explicitly and leads to an 
effective Schrodinger equation in momentum space, 

# eff |*o)=£ SP |*o>, (18) 

H cS = £ e ^<Aa + %£(# T £*j + h.c.) , (19) 



with 



k.a 



4 lHi ' • (2°) 

id 

The one-particle Schrodinger equation (|18H is readily solved, 

H cS = (h+ h kfl + ht tl h ktl ) + const. , (21) 

fe 

by means of a Bogoliubov transformation 

Cfc.T = u khk,a + v kh~t k<1 , 

ct k l = -v k hk,o + Ukht kt i ■ (22) 
Here, the real coefficients u k , Vk and the energies Ek are determined by the equations 
E k = ek(u 2 k ~ ujfe) - 2r/ s u k v k , (23) 

and 

2e k u k v k +%(u 2 k -vt)=0 , ul + v 2 k = l, (24) 
which are solved by 



E k = S ign(e k )^4 + ^ = sign(e°)^(Qe« + 7? n ) 2 + ^ , (25) 



i r. Th 
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The one-particle state |\&o) is chosen as the ground state of (|21[l 

i*o>= n ( 2g ) 

k(E k <0) 

In it was shown that the eigenvalues (|25|l in the effective Hamiltonian (|2f (I can be 

interpreted as quasi-particle energies that might be measured in ARPES experiments. 
Therefore, the parameter r] s describes the superconducting gap and Q is a measure 
for the band-width renormalization. 

When we introduce the uncorrelated density of states 

0® = (27) 

k 

the elements of the local density matrix (JSJ can be written as 

A = | / dW{l) 1 == , (29) 

and the one-particle energy reads 



4. Results for half band-filling 

We assume that the density of states is symmetric and focus on the half-filled case in 
which certain simplifications occur. By setting n n = we ensure that 2no = n = 1 
and A = (c%c^\)g = Ao. The remaining numerical task is the minimisation of the 
variational energy 

E vai = Q(A ,X cl )E° iri (A ,r )s ) + Ud(A ,X d ) (31) 

with respect to n s and A^ where Ao is given by equation (|29|l . In a pure mean-field 
BCS theory we have to set Ad = 1 and minimise the energy only with respect to rj s . 
Note that for the numerical minimisation we found it more convenient to work with 
the two variational parameters d and y = r/ s /(2Q). 



4-1- Comparison with perturbation theory 

For small values of U, and, correspondingly, small values of i] s and Ao, the 
minimisation can be carried out analytically. To leading order in rj s and U the 
variational energy reads 

n 2 



E V M « e - M - D(0)^ln ) + Fn, /)(())-,/; 



ln(^ 



(32) 

with the band-width and the bare kinetic energy 

e a = 2 [ eD(e)de . (33) 
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The difference between a Gutzwiller and a BCS wave function shows up solely in the 
respective values of the coefficient ol\j in (|3*2"|) . It is au = 1 in standard BCS theory and 
ajj = 1 — (3U)/(l6e ) < 1 for a Gutzwiller-correlated BCS state. The minimisation 
with respect to ?7 S then yields 

W ( 1 1 \ 

% = T exp (-2- W^J- (34) 

In the limit f7 — > the order parameter A ~ ?7 S in the Gutzwiller theory is therefore 
just renormalized in comparison with the respective BCS parameter 

A GW /A BCS = r = exp (__^ (35) 

An exact calculation of A by means of perturbation theory ^| E] yields the same 
functional form, eq. (|35|1 . however, with a different parameter C. It is C = 3/(16|eo|) 
in the Gutzwiller theory and 

O— < = 2j™ J^dede'^^l (36) 

in the exact evaluation. In the rest of our work we employ a semi-elliptic density of 
states of width W — At, 

D{e) = ±y/T^, (37) 
Z7T 

which is realized in a Bethe lattice with infinite coordination number. For this semi- 



elliptic density of states we find C somi = 9tt/128 w 0.221 and C^' = 4/(3tt) « 0.424. 
We then have 



/ 9vr 2 \ 

r scmi = exp f - — j w 0.4996 , (38) 

in comparison with the exact value r™^' sa 0.2636. 

Another quantity of interest is the condensation energy, i.e., the energy difference 
between the superconducting and the normal ground state. Both in our approach and 
in perturbation theory, one finds to leading order in U 

E CODd ^-±D(0)U 2 A 2 =r 2 EZ d - (39) 

For the semi-elliptic density of states the condensation energy is therefore overes- 
timated by a factor of ( r / r cxact ) 2 which shows that BCS theory overestimates the 
condensation energy by a factor of about 14 for a semi-elliptic density of states whereas 
Gutzwiller theory is only off by a factor of four. 

In the large U limit, the BCS energy shows an unphysical divergence of the 
condensation energy, Ef^ ~ U. Our Gutzwiller approach yields the correct \jU 
behaviour, 

£cond = ~Jj- (40) 

However, the coefficient 

r° 

5=2 / D(e)e 2 de (41) 

J — oo 

is twice as large as the exact result. This deficit is due to the poor description of the 
paramagnetic insulator in the Gutzwiller theory. 
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Figure 1. Rcnormalisation factor Q (solid line) and ratio of double occupancies 
d-BCs/d (dashed line) as a function of |f/|/t. Inset: double occupancies d (solid 
line) and decs (dashed line). 



4-2. Numerical results 

As seen in eq. I|25[) the Gutzwiller-Bogoliubov quasi-particles have a dispersion relation 
\Ek\ = J (Qe'fc) 2 + ??s )■ The band- width renormalization factor Q in the Gutzwiller 
theory deviates only slightly from the uncorrelated BCS value Q = 1, as shown in 
Fig- K has a minimum q m i n ~ 0.95 for medium-size correlation strength \U\ ~ 3. 
In the limit U — ► — oo the BCS wave function itself has the optimum expectation 
value d ~ 1/2 of double occupancies such that the Gutzwiller correlator P cannot 
achieve any further energy gain. Therefore, the renormalization approaches unity in 
this limit. In Fig.^we also plot the average double occupancy in Gutzwiller and BCS 
theory (see inset) as well as the ratio of these two quantities. For small values of \U\ 
the double occupancies show a qualitatively different behaviour in both approaches. 
In Gutzwiller theory we find the well-known linearity (d — 1/4) ~ U whereas BCS 
theory gives the typical exponential dependence, (c^bcs — 1/4) ~ Aq ~ exp[— (1/|J7|)]. 
The BCS value c?bcs is smaller than d only for small values of \U\. For larger \U\ the 
limited flexibility of the BCS wave function leads to an overshooting of the double 
occupancy. 

In Fig. we show the values of the superconducting order parameter A in both 
theories as a function of U. The inset of this figure shows the ratio of the two 
superconducting gaps. Note that in BCS theory we have rjf cs = |C/|A BCS . The 
order parameters are quite similar in both approaches apart from small values of |J7| 
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Figure 2. Superconducting order parameters for the Gutzwiller wave function 
(solid line) and the BCS wave function (dashed line) as a function of |C/|/t. Inset: 
Ratio of the superconducting gaps in the Gutzwiller and the BCS theory. 



where 



lira 



= lim 



iBCS 



(42) 



|t/|— o A BCS \u\->o 77^ 

Other than the order parameter, the superconducting gaps in both theories differ 
significantly over a large range of correlation parameters \U\. Equation (|35l) shows 
that BCS theory fails when the density of states at the Fermi energy -D(O) becomes 
small. 



5. Summary 

As a first application of our recently developed Gutzwiller theory for superconducting 
systems |Hj we investigated the attractive Hubbard model. Our variational approach 
leads to the diagonalisation of an effective one-particle Hamiltonian, very similar 
to the corresponding mean-field BCS theory. Quantitatively, however, the resulting 
quasi-particle band structure differs remarkably from the BCS band energies. This 
observation holds, in particular, for the superconducting gap in both theories. 

In the limit of small Coulomb interaction \U\/t, a comparison of our approach with 
perturbation theory shows that Gutzwiller theory reproduces the renormalization of 
the superconducting order parameter qualitatively. Quantitatively, the results depend 
on the bare density of the states of the system. For the semi-elliptic density of states, 
for example, our renormalization of the gap is too small by a factor of about two. 
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The condensation energy, i.e., the energy difference between the paramagnetic 
and the superconducting ground state, is another quantity which is notoriously 
overestimated in BCS theory This holds, in particular, for large U where the 
condensation energy should approach zero but diverges in BCS theory. In contrast, 
our Gutzwiller theory gives the correct order of magnitude for the condensation energy 
for all Coulomb strengths. 

Apparently, the investigation of superconducting systems with Gutzwiller wave 
functions leads to a significant improvement over on simple BCS-type calculations. 
Therefore, our approach should be a useful tool for the investigation of more realistic 
model systems that exhibit superconductivity. 
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